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ABSTRACT 

We report on additional observations of degree scale anisotropy at 90 GHz 
from the Amundsen-Scott South Pole Station in Antarctica. Observations 
during the first season with the Python instrument yielded a statistically 
significant sky signal with an amplitude of AT/T ~ 3.5 x 10~ 5 for a Gaussian 
autocorrelation function model with a coherence angle 9 C = 1°. In this paper 
we report the confirmation of that signal with data taken in the second year, 
and on results from an interleaving set of fields. Using the entire data set, we 
find AT/T = v^Co = 2-8-1'} x 10~ 5 for the Gaussian autocorrelation model 
mentioned above, and AT/T = yje e (£ e + 1)Q/(2tt) = 2.1±° 7 5 x KT 5 for a 
band-power estimate, where £ e = 93 is the effective center of our window 
function. The stated errors represent a 68% confidence interval in the likelihood 
added in quadrature with a 20% calibration uncertainty. 

Subject headings: cosmic microwave background — cosmology: observations 
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1. Introduction 



Anisotropies in the Cosmic Microwave Background Radiation (CMBR) contain a 
wealth of information about the conditions and processes that led to the formation of large 
scale structures in the universe. Whether the correct model of the evolution of the universe 
involves baryons, exotic cold dark matter, massive neutrinos, topological defects, a standard 
ionization history or reionization at a modest redshift, the content and history of the 
universe affect the observed anisotropies of the CMBR (see [White, Scott and Silk (1994)1 for 



a review). Constrained by the measured level of anisotropy on large angular scales ( |Bennett 



et al. 1994| , |Canga et al. 1993]) , different models predict varying levels of anisotropy on 



smaller scales. Several observing teams ( |Cheng et al. 1994] , |de Bernardis et al. 1994| , pevlin| 



et al. 1994] , |Cundersen et al. 1995| , [Netterfield et al. 1995| ) including our own have reported 



detections of anisotropy on angular scales near 1°. Improvements in these measurements 
should help discriminate between the current models of large scale structure formation. 

The Python instrument is designed to search for anisotropies at an angular scale 
near 1°. The first observations (hereafter Pyl) with Python were made between 1 and 15 
January 1993 at the Amundsen-Scott South Pole station ( pragovan et al. 1994| , hereafter 



Paper I). Statistically significant signals were detected in observations at high Galactic 
latitude (|6| > 49°). 

We report here the results of a second set of observations (hereafter Pyll), made 
between 12 and 23 December 1993 from the same site. The bulk of the observing time 
in the second season was spent on the Pyl field ("Field A"), in an effort to repeat that 
measurement. The remainder was spent on a set of interleaving spots ("Field B"). The 
observation of the Field A spots serves as a further check for the effects of interfering signals 
that would differ from one year to the next, including effects from atmospheric emission, 
cosmic ray hits in the bolometers, and radiofrequency interference. Additionally, the Sun 
was approximately 20° further away from the observed fields during the Pyll season than it 
was for Pyl. 



2. Instrument 



The Python instrument (Paper I, Ruhl 1993 ) consists of a 0.75 m diameter off-axis 
parabolic telescope that couples radiation into a 2 x 2 array of bolometric detectors. The 
beam response is well approximated by a Gaussian with a full width at half maximum 
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of (0.75 ± 0.05)°. Radiation from the sky is first reflected off a vertical flat before being 
focused into the cryostat by the primary. Rotation of the flat about a vertical axis 
moves the detector response horizontally across the sky. The telescope is mounted on an 
azimuth-elevation mount and is surrounded by a large shield that protects the instrument 
from being illuminated by the Sun or Earth. 

Inside the cryostat there are four corrugated feed horns at the focal plane of the 
primary. Radiation that enters a feed horn passes through a set of single-mode waveguide 
filters that define the passband at A = 3.3 mm before reaching the bolometric detectors. The 
detectors use a layer of bismuth evaporated on a diamond wafer as the absorbing element, 
and a chip of neutron transmutation doped germanium as the thermistor element. They are 
operated at 50 mK, cooled by a 3 He guarded adiabatic demagnetization refrigerator ( |Ruhl 



|fc Dragovan 1992] ). An additional bolometer (the "dark channel"), mounted on the cold 



stage but kept in a sealed cavity, acts as a monitor for extraneous pickup. 

New bolometers, constructed using the same methods as the originals, were installed 
for the Pyll observing season. However, two of the four new optical channels did not 
work well. One channel was a factor of 3 less sensitive than typical (6 mK-^/s instead of 
2 ulK-a/s); the other was unusable due to an electrical problem in the dewar. Fortunately, 
one good channel was on the upper row, and the other was on the lower, making possible 
the confirmation of the bulk of the Pyl data. 

The two most significant changes to the telescope between Pyl and Pyll were the 
improved balancing of the chopping mechanism for the external switching flat, and the 
installation of a microwave-absorbing guard ring around the primary. 



3. Observations 



The Sun, Moon, and two sources in the Carinae nebula were used as absolute pointing 
references each season; additionally, Venus was used during the Pyl season. We estimate 
our absolute pointing accuracy to be ±0.1° each year, or ±0.15° for the relative pointing 
accuracy of the Pyl and Pyll datasets. 

The same observing strategy was used for Pyl and Pyll, and is described in Paper I. 
The combination of a fast 3-beam chop (2.5 Hz full cycle) and slower beamswitching (10 to 
30 seconds per position) yields a 4-beam response on the sky, 

AT " = ~\ T i + ~ ~ A T J+2 + l -T 3+3 , (1) 
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where the T^-'s are the antenna temperatures of patches on the sky separated by 2.75° along 
a horizontal line. A scan consists of measuring ATj° three times successively (k = 1, 2, 3) at 
each of 7 positions j on the sky. The left and right hand channels in a given row of the 
array measure many of the same 4-beam patterns; the spots measured by the left hand side 
are given by j = 1, 7, while those made by the right correspond to j = 2, 8. 

The time-ordering of the right and left-hand 3-beam measurements ("stares") are 
reversed from one value of j to the next. For the first 4-beam pattern (j = 1 for the left 
channels, j = 2 for the right), the left stare is measured first. For the second 4-beam 
pattern the right stare is measured first, and so on. This causes a drifting 3-beam offset to 
appear as an oscillating 4-beam signal on the sky. However, the measured linear drifts in 
the 3-beam offset within each scan are both small and symmetric about zero, making this 
effect unimportant. The average 3-beam offsets are less than 1 mK for all the Pyl and Pyll 
channels, and the average slopes are within 2a of zero, with an error of a ~ 1/iK/stare. 



4. Calibration 



The Pyl observations were calibrated by a combination of elevation scans and 
the placement of known-temperature blackbody and low emissivity foam loads in the 
optical path, as described in Paper I and |Ruhl et al. 1995| . The Pyll observations were 
calibrated using the foam loads; the loads were in turn calibrated by placing a large 
liquid-nitrogen-cooled load beneath the dewar and alternately switching the foam and two 
known-temperature loads into the optical path. This calibration of the foam loads was 
done in situ on the telescope during the Pyll season, rather than in a 300 K room, reducing 
possible systematic effects in the measurement. The DC gain measured using the foam 
loads was converted to an AC gain by comparing the signals seen while switching on and 
off a foam load slowly (0.1 Hz) with those seen while switching at the frequency used for 
the observations (5 Hz). 

The relative calibrations of the various bolometer channels (four channels for Pyl, 
three for Pyll) were checked using two sources in the Carinae nebula; the results from all 
seven channels lie within 15% of the average. The two sources in the Carinae nebula are at 
(a, 5) = (10 h 44 m , -59.64°) and (10 h 33 m , -57.95°) in J1994 coordinates. The mean signals 
from them are (9.0 ± 0.3) and (4.2 ± 0.1) mK in our beam, respectively, with the errors 
representing only the statistical uncertainty. The first source was also used as a pointing 
reference. 
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The statistical error on the determination of the average gain is 5%. Our estimates of 
possible systematic effects are at the level of < 20%. We therefore adopt a gain uncertainty 
of 20%. 



5. Analysis 

The data analysis presented here differs slightly from that used for the previously 
reported Pyl results; using the new method does not significantly change those results. We 
describe here the new method, which is used to arrive at all the results in this paper. 

The output of each detector is sampled at 100 Hz and demodulated in software, 
synchronously with the motion of the chopping flat. Two lockin demodulations are used; 
one (the "optical phase") is maximally sensitive to the 3-beam signal from the sky, while 
the second (the "quadrature phase") is shifted by 90° from the optical phase and should 
have no response to a stationary sky signal. We confine our discussion to the optical phase 
data, giving results from the quadrature phase where relevant. 

Within each of the 42 stares in a scan, the mean and variance erf of the lockin values 
are calculated. The means are then combined in a pairwise fashion into 4-beam values AT* 
as described above, and an average variance a% is calculated from the 42 values of a\.. Thus 
erf is a measure of the noise which contains drifts only within 10 to 30 second long stares. 

As previously described, each channel measures a 4-beam value AT* three times in 
succession (k = 1,2,3) for seven sky positions j. From these we find a mean for each 
channel at each position j, and an average error on those means, a m . Drifts up to the time 
separating observations of successive 4-beam positions (roughly 3 minutes), are included in 
a m . The variance from a celestial signal does not contribute to a m , and a m can therefore 
be used as an unbiased statistic for cutting scans. 

The noise estimates described above are calculated for each channel in every complete 
scan. We use as and a m , normalized to a "stare" duration of 30 seconds, on a channel 
by channel basis to remove scans which have been contaminated by excessive noise. The 
first cut removes scans with as greater than 1.5 times the peak of its distribution. This 
procedure removes 20 to 25% of the scans for each channel in Pyl, and 5 to 10% of those in 
Pyll. The second cut removes scans with high values of a m ; the value at which scans are cut 
is placed so as to minimize the average final errorbar in the binned data. This step removes 
19 to 27% of the scans for Pyl, and 4 to 5% of those in Pyll. The cuts remove more Pyl 
data than Pyll because of the poorer weather during that season. In all, approximately 50 
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hours (depending on the channel) of Pyl data passed the cuts. For Pyll, roughly 24 hours 
of data on the overlapping field (Field A), and 9 hours on the interleaving field (Field B) 
remained after the cuts. 

The average value of as from all scans that pass the cuts is a good indicator of chopped 
detector noise. For Pyl, the noise in the four channels is 2.5, 2.9, 2.1 and 1.9 niK -y/s; 
For the three operational channels of Pyll, it is 2.1, 1.3 and 5.7 mK -y/s. These values 
are in good agreement with those found from the quadrature values of both as and a m . 
The optical phase value of a m is higher (for the sensitive channels), due to atmospheric 
contamination. For Pyl it is 4.1, 4.2, 4.0 and 3.9 mK -y/s; for Pyll these values are 2.9, 2.7 
and 5.6 mK -y/s. 

For sky positions observed by more than one channel, the 4-beam averages from 
neighboring channels are combined into a single value within a scan by forming a weighted 
average of the left and right-hand channels. The average value of over all uncut scans 
is used as the weight. If one of the channels is cut by the previously described procedures, 
its neighbor is cut as well. 

After this treatment, a statistical mean and error are calculated from the uncut scans 
at each of the sky positions j. The errors are within 10% of those expected given the 
distribution of a m , indicating there is little if any additional atmospheric contamination. 
These means and errors, after multiplying by 1.07 to correct for the estimated atmospheric 
absorbtion in our band, are the final values and errors for the 4-beam temperature 
differences on the sky. 



6. Results 

The results of the analysis of the Pyl data and of the overlapping portion (Field A) of 
the Pyll data set appear in Figure |I[ The agreement is good; taking the difference of the 
two data sets [ (Pyl — PyII)/2 ] gives a result that agrees well with zero signal (reduced 
X 2 = xt = 12./15). The weighted mean of the two data sets is shown as the filled circles 
in Figure ^ here the agreement with zero signal is very poor (xl = 191./16). The contrast 
between the x 2 ' s for the summed and differenced data sets indicates that the signal in the 
two data sets is the same. 

A set of patches (Field B) that interleaves those from Field A was also observed during 
the Pyll season. The locations of the Field B beams are found by moving one half of a 
chopper throw in negative right ascension (1.38° on the sky) from the beams that make up 
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Field A. Analyzing the Field B data leads to the values shown as open squares in Figure 0. 
Less time was spent observing this field, and the errors are larger than on Field A. The 
agreement with zero sky signal is good for Field B (xf, = 14. /15). 

The results shown in Figure || were measured in a single frequency band, so we cannot 
spectrally discriminate between CMBR anisotropics and the various possible foregrounds. 
In Paper I we discussed the expected levels of galactic foreground contamination, which 
are smaller than the signals seen. Proceeding under the assumption that these signals 
represent fluctuations in the CMBR, we multiply the values plotted in Figure |] by 1.24 to 
convert them from Rayleigh- Jeans temperature units to temperature differences in a 2.73 K 
blackbody, and use the integrated likelihood function to form confidence interval estimates 
for the sky signal in two models with Gaussian CMBR fluctuations. For both analyses 
we use the full correlation matrix including off-diagonal elements describing correlations 
present because of the theoretical model and those induced by the observing strategy. 

The first model consists of a sky with a Gaussian autocorrelation function (GACF) 
with a coherence angle of 9 C = 1°. We set limits on \/Cq, the rms amplitude of the 
fluctuations in AT/T. The uncertainty in beamwidth leads to less than a 2% error in \/Co. 
The GACF limits can be converted to an estimate of the flat band power sampled by our 



window function. Following the recipe of [Bond 1994] , we find for our window function 



AT Band p ower = 0. 73 AT G acf- The effective center of our window function for a flat band 
power spectrum lies at £ e = 93. 

We also use a model with an uncorellated sky, C(9) = 5(1 — cos 9), to find the rms of 
our data set, from which we derive another band power estimate. Correlations due to beam 
overlap are less than 1%, and we ignore them. This analysis gives results that are within 
a few percent of those found by converting the GACF limits, indicating that the flat band 
power is insensitive to the form of the corellation function used in the likelihood calculation. 

The results for the various combinations of data sets and analyses are given in Table |I[ 
The confidence intervals quoted in the table do not include the 20% calibration error. The 
first four entries in Table [I] show that the signal seen in the first season appears again in the 
second. The last entry in Table [I] gives our final result using the data from both seasons 
and both fields. 



7. Conclusions 
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After detecting a signal in our first year of observations, we re-observed the same 
portion of the sky and detected the same signal. This provides further evidence that the 
signal is celestial rather than systematic in nature. 

Additionally, an interleaving set of spots was observed during the second season. 
Likelihood analyses on the entire data set are used to derive two estimates of the sky signal. 
First, we find 79l^g/xK for the total sky rms for a Gaussian autocorrelation model with a 
coherence angle of 1°. Second, we find a band power of 57l?VK, using an autocorrelation 
fuction given by C{9) = 5(1 — cos#). The stated errors are 1 a limits that include the 20% 
calibration uncertainty added in quadrature with the likelihood-derived errors. 

We thank the Antarctic Support Associates staff at the South Pole for making a 
successful season possible, and Ted Griffith, Bob Pernic, and Bill Vinje for valuable 
assistance there. This research was supported by the James S. McDonnell Foundation, 
PYI grant NSF AST 90-57089, and the National Science Foundation under a cooperative 
agreement with the Center for Astrophysical Research in Antarctica (CARA), grant NSF 
OPP 89-20223. CARA is an NSF Science and Technology Center. JR was supported by the 
McCormick Fellowship at the University of Chicago. 
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Dataset Field Band Power GACF 



4 = 93 e c = i° 



Pyl 


A 


66±?I 


90^ 


Pyll 


A 


68+?? 




(Pyl&Pyll) 


A 




93±?| 


(PyI-PyII)/2. 


A 


(3,19) 


(5,26) 


Pyll 


B 


(16,59) 


(23,81) 


(Pyl&Pyll) 


A&B 


571? 





Table 1: Results of likelihood analyses. 

Values are in units of /iK. Detections are quoted at the maximum in 
the likelihood, with 16% of the integrated likelihood above the upper 
errorbar, and 16% of the integrated likelihood below the lower one. For 
datasets with no significant detection, the value which maximizes the 
likelihood is given, along with a 95% upper limit from the integrated 
likelihood. The 20% calibration uncertainty is not included in these 
errors. The signals from the dark channel and from the quadrature 
phase are consistent with zero. 
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Fig. 1. — Results from Pyl and Pyll for Field A. The Pyl results are plotted as the solid 
circles, and the Pyll results are shown as the open squares. They are offset horizontally 
for clarity. The 4-beam response pattern is shown in the top panel, aligned with the point 
plotted at a = —12°. The temperature axis is plotted in Rayleigh- Jeans units; to convert to 
thermodynamic temperature differences for a 2.73 K blackbody, multiply by 1.24. 



Fig. 2. — Final results for both fields. The weighted mean of Pyl and Pyll for Field A are 
shown as the solid circles. The open squares show the results of the observations of Field B 
done during the Pyll season. The temperature axis is plotted in Rayleigh- Jeans units; to 
convert to thermodynamic temperature differences for a 2.73 K blackbody, multiply by 1.24. 
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